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137  MHz.  Equations  are  developed  which  yield  scintillations  at  this  fre- 
quency as  a function  of  local  time,  magnetic  index,  solar  flux,  and  month  of 
the  year.  The  conceprt  is  to  predict,  at  sub-auroral  and  auroral  latitudes, 
scintillation  at  this  frequency  and  at  higher  frequencies.  The  aim  is  to  give 
to  users  an  indication  of  expected  scintillation  excursion  when  one  predicts 
magnetic  activity  and  solar  radio  flux.  The  equations  developed  have  been 
checked  with  one  additional  set  of  observations  and  are  being  checked  with 
additional  data. 

Phase  2 of  this  model  will  incorporate  geometrical  terms  to  take  account 
of  the  propagation  angle  of  the  observer  vis-a-vis  the  irregularities  and  will 
allow  for  frequency  dependence  to  be  ascertained.  In  addition,  it  is  expected 
in  Phase  2 to  validate  the  model  with  additional  sets  of  data,  extending  the 
model  to  auroral  latitudes  greater  than  63®ipind  to  polar  latitudes. 
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A High-Latitude  Empirical  Model 
of  Scintillation  Excursions:  Phase  1 

I 

1.  INTRODUCTION  i 


In  order  to  provide  a realistic  series  of  values  to  determine  scintillation  oc- 
currence at  sub -auroral  and  auroral  latitudes,  the  Trans -Ionospheric  Branch  of  [ 

the  Space  Physics  Laboratory  in  collaboration  with  Loglcon,  Inc.  initiated  a model  ! 

development.  j 

The  data  base  was  to  be  a series  of  measurements  of  the  scintillations  of 
synchronous  satellite  beacons  at  137  MHz  extending  over  several  years.  The  initial 
model  was  to  provide  mean  scintillation  excursions  as  a function  of  time  of  day, 
month,  solar  flux,  and  magnetic  index.  With  forecasting  of  solar  flux  and  magnetic 
index,  a user  in  the  latitude  range  covered  would  have  an  indication  of  mean  scintilla- 
tions to  be  expected. 

The  information  would  be  useful  in  estimating  scintillations  at  higher  frequen- 
cies such  as  250  to  400  MHz  and  at  other  geometries  in  the  sub -auroral  and  auroral 
ovals.  In  addition  the  data  could  be  used  to  forecast  clutter  problems  on  HF  back- 
scatter  systems  that  show  difficulties  from  the  same  F layer  irregularities  that 
produce  scintillations. 


(Received  for  publication  15  September  197  6) 
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2.  THK  DATA  BASE 


Observations  of  the  scintillations  of  the  ATS -3  satellite  beacon  at  137  MHz 
form  the  data  base.  The  data  was  reduced  by  the  method  outlined  by  Whitney 
et  al.  Scintillation  excursions  in  dB  for  15 -min  periods  are  used;  these  approxi- 
mately correspond  to  the  1%  percentile  levels  in  the  cvimulative  amplitude  probabil 
ity  distribution  function  for  the  15-min  period.  The  dB  excursions  are  convertible 
into  indices  by  an  empirical  relationship.  The  three  stations  are  situated  near 

the  70°  West  meridian,  and  their  propagation  paths  to  the  ATS -3  satellite  are  shown  in 

FigUi  e 1.  Their  geomagnetic  and  geometrical  parameters  are  given  in  Table  1. 


Thule 


'narssarssuao 


Figure  1.  Observing  Stations  (ATS -3 
at  70°W,  LES-6  at  39°W,  ATS-5  at 
105°W).  Thule  data  was  not  used  in 
Phase  1.  X denotes  the  350  km 
intersection  point  to  the  satellite. 
Propagation  paths  to  ATS -3  are  shown 
as  well  as  those  to  LES-6  and  ATS-5. 
Only  data  from  ATS-3  was  used  in 
Phase  1.  Other  data  will  be  used  in 
Phase  2 


GOOSE  BAY 


SAGAMORE  HILL 


50*  INVARIANT 
LATITUDE 


Table  1.  Geomagnetic  and  Geometrical  Parameters  for  ATS-3 


Inv. 

Lat. 

El. 

Az. 

lono. 

Z.  A. 

63.2° 

18.0° 

208° 

64° 

60.3° 

28.  8° 

191° 

56° 

53.  5° 

40.  9° 

178° 

46° 

3.  DATA  BASE  CREATION  AND  MODELING 


Fifteen -minute  samples  of  scintillation  index  (dB)  data  for  the  ATS -3  satellite 
from  the  above  three  stations  have  been  augmented  into  complete  data  bases  on  tape. 
Data  associated  with  each  15 -min  sample  include:  Station  ID.  Satellite  ID,  Frequen- 
cy (137  MHz),  Date  and  UT,  Scintillation  Index  (dB),  3 -hourly  magnetic  index  (K  ). 
2695  MHz  solar  flux,  (S^.),  4995  MHz  solar  flux,  sub -ionospheric  latitude  and  longi- 
tude. and  corresponding  geomagnetic  latitude  and  magnetic  local  time. 

Narssarssuaq  9/17/68  - 9/1/74  146, 700  + samples 

(3oose  Bay  1/1/72  - 12/31/74  7 1. 000  + samples 

Sagamore  Hill  12/1/69  - 11/30/74  148, 000  + samples 

Ctoose  Bay  data  for  1974  cover  mainly  November  and  December. 

4.  THE  FORCING  FL’NCTIONS 


In  the  initial  development,  termed  Phase  1,  a relatively  simple  equation  was 
developed  to  reproduce  the  mean  scintillation  excursions  at  each  station. 

The  forcing  parameters  included  the  following: 

(a)  Planetary  magnetic  index,  Kp.  It  has  been  shown  that  at  high  latitudes 

(which  in  this  case  include  the  Narssarssuaq  and  (3oose  Bay  observations  and.  dur- 
ing severe  magnetic  storms,  the  Sagamore  Hill  observations)  K is  a forcine  func- 
*.2  P ® 

tion. 

O 

(b)  Early  work  by  Aarons  et  al  showeJ  that  solar  flux,  even  when  divorced 
from  magnetic  index  variations,  had  an  effect  on  scintillation  behavior.  Utilizing  the 
2695  MHz  solar  measurements  at  Sagamore  Hill,'^  we  separated  observations  into 
three  regimes  of  solar  flux  units  (S^,);  (1)  0-95,  (2)  96-120,  and  (3)  121  and 
greater.  As  will  be  shown  scintillations  do  not  in  every  month  increase  with  increas- 
ing solar  flux  but  vary  as  a function  of  season. 

(c)  The  seasonal  parameter.  A very  dramatic  minimizing  of  diurnal 
effects  is  shown  in  the  Narssarssuaq  observations  over  the  winter.  ^ A similar 
pattern  is  seen  on  the  (3oose  Bay  records.  This  would  have  to  be  reflected  in  the 
model. 

2.  Aarons,  J.  and  Allen,  R.S.  (197  I)  Scintillation  boundaries  during  quiet  and 

disturbed  magnetic  conditions,  J.  Geophys.  Res. . 17 0-177. 

3.  Aarons,  J.,  Whitney,  H.  E. , and  Allen.  R.S.  (1971)  aobal  morphology  of 

ionospheric  scintillations,  Proc.  IEEE. . ^:159. 

4.  Solar -Cteophysical  Data.  lER-FB  289-304,  305-364  Part  I,  U.  S.  Department  of 

Commerce  (Boulder,  Colorado,  U.  S.A.  80302). 

5.  Basu,  Sunanda  (1975)  Universal  time  seasonal  variations  of  auroral  zone  mag- 

netic activity  and  VHF  scintillations,  J.  Geophys.  Res.,  80:4725-4728. 
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(d)  The  diurnal  pattern  would  have  to  be  shown.  This  would  include  the 
variations  as  a function  of  local  time  under  magnetically  quiet  and  disturbed  condi- 
tions as  modulated  by  the  seasonal  function. 

Analyses  were  conducted  separately  for  each  station.  The  data  were  further 
partitioned  into  12  months,  7 K^,  3 Solar  Flux  (2.7  GHzi,  and  24  UT  ranges.  The 
date,  Kp,  S^.  and  SI  readings  were  averaged  in  each  block.  A compact  file  was 
thus  made  available  for  high  speed  iterative  modeling  studies.  The  seven  Kp 
ranges  are  0-1,  1+  to  2,  2+  to  3,  3+  to  4,  4+  to  5,  5+  to  6,  6+  and  up.  The  three 
Sj,  ranges  are  0 to  95,  96  to  120,  121  and  up.  Tables  of  the  averaged  SI  are  pro- 
vided for  each  of  the  stations. 

Out  of  a maximum  possible  6048  blocks  (12  x 7 x 3 x 24),  the  averaged  files 
comprise  the  following: 

Narssarssuaq  4985  blocks 

Goose  Bay  4217  blocks 

Sagamore  Hill  5065  blocks. 

The  empty  blocks  generally  correspond  to  the  highest  two  Kp  ranges,  that  is,  5+ 
and  up,  and  occasionally  the  highest  S^  range. 


5.  MODKLING 

An  extensive  search  was  conducted  to  derive  empirical  models  of  SI  for  each 
of  the  three  stations  (Figure  2a).  Analytical  forms  of  the  model  were  preferred  to 
ensure  smooth  transitions  as  a function  of  the  known  driving  functions,  namely, 
day  of  year,  Kp,  .S^  and  universal  time.  These  forms  also  permitted  use  of  re- 
gression techniques  for  least  squares  fitting  tc5  tiTC"averaged  data  file.  In  the  course 
of  the  search  for  improved  fits,  special  characteristics  of  the  data  were  noted  which 
suggested  elaborations  of  the  model  form.  Examples  are  the  delayed  peak  in  the 
diurnal  SI  variation  with  higher  Kp,  the  seasonal  effect  on  diurnal  variation  ampli- 
tude relative  to  the  average  SI,  the  seasonal  effect  on  influence  of  Kp  and  S^,  and 
the  need  for  higher  harmonics  to  represent  the  diurnal  variation. 

The  models  for  average  SI  (dB)  for  the  three  stations  are  provided.  In  the 
model  equations,  the  cosine  arguments  (in  radians)  that  include  the  terms  for  day 
number  (DA)  and  hour  (HR)  assume  a multiplicative  term  of  2ir/365  and  2ir/24  re- 
spectively. Though  not  shown  in  the  equation  for  reasons  of  convenience,  these 
must  be  included  in  the  calculations. 

Figure  2b  presents  the  model  for  each  station  in  terms  of  local  time.  The  dif- 
ference between  universal  time  and  local  time  at  the  sub-ionospheric  point  (350  km) 
was  taken  to  be  3.4  hr  at  Narssarssuaq,  4.  1 hr  at  Goose  Bay,  and  4.7  hr  at  Saga- 
more Hill.  Subsequent  work  directed  towards  developing  bridging  models  used  local 
time  rather  than  universal  time  as  the  argument. 
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IN  SOLAR  FLUX  UNITS 


In  Figures  3 to  6 typical  comprehensive  comparative  plots  of  the  averaged 
data  and  the  model  are  provided  for  each  station.  The  model  predictions  used 
the  actual  averaged  data,  and  Sj.  for  each  hour,  and  are  therefore  absent  when 
data  are  absent.  These  best  fit  models  may  occasionally  predict  small  negative 

Si's;  these  are  made  to  asymptotically  approach  zero  by  replacing  a value  V which 
2 (V - 85 ) 

is  less  than  0.  5 by  e ’ . This  ensures  continuity  in  the  predictions  near  0.  5 

and  introduces  minimal  distortion  since  only  the  very  low  scintillations  are 
adjusted. 


6.  THK  equations  AND  THE  DATA 


To  illustrate  the  behavior  of  the  data  and  of  the  best  fitting  model  equations. 
Narssarssuaq  data  is  shown  in  Figure  3 for  the  month  of  March.  One  can  note, 
for  example,  along  the  left-hand  side  for  the  solar  flux  set  of  0 to  95,  that  scin- 
tillations increase  with  increasing  Kp.  The  dotted  line  is  the  model  equation  and 
the  continuous  line  connects  1-hr  means  of  scintillation  excursions. 


NARSSARSSUAQ  MONTH  3 
96-120 


SOLAR  FLUX  0-95 


Figure  3.  Typical  Comparative  Plots  of  Averaged  Data  and  the  Model 
Narssarsuaq,  March 


Along  any  series  of  boxes  with  identical  one  can  note  that  the  model  and 
the  data  indicate  that  the  general  trend  is  for  the  scintillation  excursions  to  in- 
crease with  increasing  solar  flux. 

In  the  Narssarssuaq  data  for  October  (Figure  4)  it  can  be  noted  that  with 
increasing  K (for  a constant  solar  flux)  scintillation  excursions  increase,  and  1 
scintillation  excursions  decrease  with  increasing  solar  flux  (and  constant  K ). 


NARSSARSSUAQ  MONTH  10 
96-120 


SOLAR  FLUX  0-95 


Figure  4.  Typical  Comparative  Plots  of  Averaged  Data  and  the  Model 
Narssarssuaq,  October 


GOOSE  BAY  MONTH  3 

SOLAR  FLUX  0-95  96-120  121-UP 


Figure  5.  Typical  Comparative  Plots  of  Averaged  Data  and  the  Model; 
Goose  Bay,  March 


The  behavior  of  the  Cioose  Bay  data  (solid  line)  and  best -fitting  model  (dashed 
line)  for  March  are  shown  in  Figure  5.  Within  each  particular  solar  grouping 
scintillations  increase  with  increasing  K^.  For  the  same  grouping,  the  scintilla- 
tion excursions  increase  with  increasing  solar  flux. 

Behavior  of  the  Sagamore  Hill  data  and  model  for  March  are  shown  in  Figure 
6.  For  the  low  solar  flux  group  (0  to  95),  a very  slight  increase  in  scintillations 
associated  with  increasing  is  seen.  The  same  is  true  of  the  solar  flux  group 
96  to  120  while  for  high  solar  flux,  the  increase  with  becomes  more  noticeable. 

In  general,  within  a particular  grouping,  scintillation  index  increases  with 
increasing  solar  flux. 

Although  equations  were  developed  for  Sagamore  Hill,  Goose  Bay,  and  Nars- 
sarssuaq  observations,  it  is  not  certain  whether  or  not  to  join  the  Sagamore  Hill 
observations  to  the  high  latitude  data  in  the  model  which  will  be  finally  developed. 

The  sub-auroral  to  middle  latitudes  traversed  (at  53°  invariant  latitude)  present  a 
different  morphology  from  the  high  latitude  ionosphere  except  during  periods  of 
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Figure  6.  Typical  Comparative  Plots  of  Averaged  Data  and  the  Model 
Sagamore  Hill,  March 


intense  magnetic  activity  when  auroral  effects  at  F layer  heights  cover  this  region 
Further  work  on  this  is  exhibited  in  Appendix  B. 
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Distribution  analyses  were  conducted  for  each  of  the  three  stations.  The  effect 
of  the  solar  flux  was  not  evaluated,  but  partitioning  by  month,  and  hour  was 
retained.  Within  each  partition  the  average  SI  was  calculated,  and  the  percentage 
of  readings  falling  in  consecutive  ranges  was  determined.  Significant  features  of 
the  distributions  have  been  obtained  by  examination  of  the  tabulations,  and  are  dis- 
cussed below.  For  selected  average  SI,  the  percentage  SI  readings  found  in  the 
ranges  0-1,  1-3,  3-6,  6-9,  9-12,  and  12  and  up  are  presented  in  Table  2. 


Aarons,  .1.  (1976)  High-latitude  irregularities  during  the  magnetic  storm  of 
October  31  to  November  1,  1972,  .1.  Geophvs.  Res.,  81:661-670. 


Table  2.  Distribution  of  SI  Readings  for  Various  Average  SI 
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Percent  Occurrence 
SI  Range 

k. 

1 

Avg.  SI 

0-1 

1-3 

3-6 

6-9 

9-12 

12  up 

1 

Narssarssuaq 

50 

45 

5 

0 

0 

0 

f 

i 

1.  0 

Goose  Bay 

55 

35 

9 

1 

0 

0 

Sagamore  Hill 

70 

15 

10 

3 

2 

0 

Narssarssuaq 

20 

55 

20 

4 

1 

0 

i. 

2.  0 

Goose  Bay 

40 

40 

15 

.3 

1 

1 

> 

Sagamore  Hill 

55 

25 

10 

5 

3 

2 

Narssarssuaq 

20 

40 

25 

10 

3 

2 

3.  0 

Goose  Bay 

30 

33 

18 

10 

5 

4 

Sagamore  Hill 

50 

20 

12 

7 

5 

6 

1 

4.  0 

Narssarssuaq 

10 

38 

28 

15 

4 

5 

1 

Goose  Bay 

20 

33 

20 

15 

5 

7 

i' 

5.  0 

Narssarssuaq 

6 

33 

30 

18 

5 

8 

\ 

Goose  Bay 

12 

30 

25 

18 

6 

9 

k' 

6.  0 

any 

3 

25 

35 

20 

6 

11 

1 

8.  0 

anv 

2 

15 

25 

20 

16 

22 

10.  0 

any 

1 

6 

15 

18 

20 

40 

13.  0 

any 

1 

1 

8 

10 

20 

i 

60 
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For  each  station,  no  marked  variations  in  the  distributions  are  evident  for 
different  months,  ranges  or  time  of  day.  Thus,  given  an  average  SI  value  for 
any  time-K^  partition  for  a specific  station,  the  expectation  of  any  range  of  read- 
ings may  be  predicted.  Differences  in  distribution  patterns  among  the  three  sta- 
tions are  primarily  due  to  the  consistently  higher  SI  readings  at  Narssarssuaq  on 
the  one  hand,  and  the  extremely  low  activity  at  Sagamore  Hill  on  the  other;  0 to 
1 dB  readings  are  rarer  at  Narssarssuaq  so  that  average  SI  of  1 to  3 dB  are  found 
with  negligible  occurrence  at  the  higher  SI  readings.  At  Sagamore  Hill  0 to  1 dB 
readings  constitute  more  than  50  percent  of  the  samples  at  any  time  so  that  an  aver- 
age of  even  1 dB  implies  the  presence  of  some  high  SI  readings.  Though  the  1 to  3 dB 
averages  usually  occur  in  the  daytime  at  Narssarssuaq  and  in  the  night  sector  at 
Sagamore  Hill,  no  particular  diurnal  influence  is  evident  on  the  distributions.  Thus, 
the  instances  where  low  average  SI  are  found  in  the  nighttime  at  Narssarssuaq  are 
devoid  of  higher  SI  readings  as  in  the  daytime;  the  instances  where  average  SI  at 
Sagamore  Hill  are  over  1 dB  in  the  daytime,  higher  sporadic  SI  readings  are  clear- 
ly the  cause. 

The  above  discussion  covers  the  latitudinal  and  SI  activity  extremes.  Goose 
Bay  lies  in  between  and  characteristics  consistent  with  the  above  arguments  are 
exhibited  in  the  distributions  for  this  station.  As  seen  in  Table  2,  below  an  aver- 
age of  6 dB.  Goose  Bay  readings  tend  to  include  more  high  SI  samples  than  Nars- 
sarssuaq. 
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8.  (;komf.tric.\i.  paramktkrs 


The  development  of  the  equations  in  the  form  given  fails  to  correct  scintilla- 
tions for  the  geometry  of  observations.  Scintillations  maximize  at  low  angles  of 
elevation  and  when  the  signal  is  observed  parallel  to  the  lines  of  force  of  the  earth's 
field.  The  complexity  of  arriving  at  a physical  model  which  determines  the  elonga- 
tion of  the  irregularities  as  a function  of  latitude,  time,  and  geomagnetic  conditions 
is  such  that  only  a simplified  solution  can  be  attempted  at  this  time. 

In  a preliminary  attempt  to  give  to  the  user  of  this  Phase  1 model  a correction 
factor  for  geometry  we  have,  from  studies  at  Narssarssuaq,  Goose  Bay,  and  Saga- 
more Hill*  presented  an  elliptical  column  model  of  individual  irregularities  of 
1 km  N-S,  2 km  orthogonal  to  the  lines  of  force  and  5 km  along  the  lines  of  force. 
Assuming  a Gaussian  irregularity  spectrum,  it  was  determined  that  the  correction 
for  the  S^  value  derived  from  each  index  would  be  a factor  of  ~2.  0 greater  at  the 
azimuth  and  elevation  of  the  ATS -3  satellite  from  Narssarssuaq  than  at  zenith. 

From  Goose  Bay  the  correction  factor  using  the  same  model  would  be  ~1.  3,  and 
from  Sagamore  Hill  it  would  be  ~1.4. 

Graphs  of  these  correction  factors  have  been  developed,  these  are  given  in 
Figures  7,  8,  and  9.  It  should  be  noted  that  the  model  equations  have  not  utilized 
these  correction  factors;  the  factors  are  given  merely  to  help  the  reader.  It  is 
expected  that  in  Phase  2 of  the  model  development  a geometrical  correction  factor 
will  be  utilized. 
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Figure  7.  Graph  of 
Correction  Factors  for 
Narssarssuaq  Under 
Assumptions  of  Elliptical 
Column  Irregularity  Model 
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7.  Singleton,  D.  G,  (1975)  Private  communication, 


S60’ 


GOOSE  BAY 


Figure  8.  Graph  of  Correction  Factors  for  (ioose  Ray  tinder 
Assumption  of  Elliptical  Column  Irregularity  Model 
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Figure  9.  Graph  of  Correction  Factors  for  Sagamore  Hill  Under 
Assumption  of  Elliptical  Column  Irregularity  Model 
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In  addition  to  the  geometrical  parameter  another  factor  for  the  frequency  de- 
pendence must  be  put  into  the  model  development.  The  frequency  correction  is 
quite  complex  and  depends  dramatically  on  whether  the  scintillations  are  in  the 
strong  scattering  or  weak  scattering  regime.  For  weak  scattering  (with  dB  excur- 
sions less  than  about  8 in  our  data  reduction  scheme)  the  expected  frequency  de- 

g 

pendence  using  S , values  is  a factor  of  1.  5.  Using  fade  depth  in  dB  (below  mean 
^ 9 

level)  the  frequency  dependence  is  1.9.  These  relationships  are  shown  in  Figure 
10  and  can  be  used  for  extrapolating  137  MHz  data  to  higher  frequencies  under  con- 
ditions of  weak  scatter. 


FREQUENCY  - MHi 

Figure  10.  Variation  of  and  Fade  Depth  with  Frequency: 
Weak  Scattering 


8.  Whitney,  H.  E.  (1974)  Notes  on  the  Relationship  of  Scintillation  Index  to  Prob- 

abilitv  Distributions  and  Their  Uses  for  System  Clesicn,  AFCR  L -Tt? -7  4 - 
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9.  Whitney,  H.  E,  , Aarons,  .1.,  and  Seeman,  D.  R,  ( 197  1 ) Estimation  of  the 

Cumulative  Amplitude  Probability  Distribution  Function  of  tonospheric 
Scintillations,  Appendix  I,  AfCRL-7  1 -0535. 
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With  strong  scattering  the  frequency  dependence  term  for  will  vary  from 
zero  to  1.  5,  The  zero  frequency  dependence  occurs  when  reduced  simultaneous 
data  at  two  frequencies  both  display  Rayleigh  distributions.  Frequency  dependence 
then  may  vary  from  zero  to  the  weak  scattering  regime  value  of  1.  9 based  on 
fade  depth. 

With  a frequency  dependence  term,  a proper  model  for  the  irregularity  config- 
uration, and  a geophysically  valid  model  of  the  effect  of  the  polar  cap,  auroral  and 
sub-auroral  terms,  it  will  be  possible  to  further  develop  this  model  so  that  it  can 
service  systems  designers  and  operators  who  have  radio  transmissions  affected  by 
high  latitude  irregularities.  This  will  be  attempted  in  Phase  2. 


10.  FUTURK  MODEL  UEVELOPMEM 

Phase  2 of  the  model  development  will  incorporate  in  its  equations,  correc- 
tions for  geometry  utilizing  newly  developed  data  on  power  law  fall  off  of  irregular- 
ity spectra  as  a function  of  latitude  and  geomagnetic  conditions.  .An  attempt  will  be 
made  to  utilize  the  data  set  out  in  this  development  by  dividing  t.he  observations  into 
weak  and  strong  scatter  regimes  to  determi.ne  frequency  dependence  for  various 
scintillation  excursions. 

Phase  " also  envisages  additional  validation  studies  of  the  model.  Data  from 
new  satellites  will  be  used  to  check  the  observations  in  the  auroral  and  sub -auroral 
regions  and  to  extend  the  observations  into  the  polar  area.  The  Goose  Bay  data  base 
was  considerably  smaller  than  the  other  stations';  it  will  be  expanded  for  the  final 
model. 

The  use  of  satellites  at  low  and  high  altitudes  transmitting  at  several  frequen- 
cies will  allow  for  high  latitude  studies. 

Multi -frequency  transmissions  will  also  allow  additional  studies  of  frequency 
dependence  to  determine  how  best  to  extend  these  137  MHz  data  to  higher  frequen- 
cies. 
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Appendix  A 

Comparison  of  Model  With  Millstone  Hill  Data 


In  order  to  validate  the  model  using  an  independent  set  of  observations  at  a 
higher  frequency  it  was  determined  that  the  Millstone  Hill  data  obtained  by  observ- 
ing scintillations  at  400  MHz  tracking  Transit  satellites  would  be  used.  The  data 
base  for  the  comparison  was  that  in  the  publication  "The  Millstone  Hill  Radar  Pro- 
pagation Study:  Scientific  Results"**’ by  .1.  Evans. 

Since  Evans'  data  was  in  the  S^  format  at  400  MHz.  a conversion  using  weak 
scatter  relationships  must  first  be  made  to  dR  at  137  MHz.  This  is  plotted  in 
Figure  A 1, 

The  measured  S^  (400  MHz)  as  a function  of  invariant  latitude  is  shown  in  Fig- 
ure A2.  These  values,  in  three  different  groupings,  at  the  invariant  latitude  of 
Narssarssuaq,  Goose  Bay,  and  .Sagamore  Hill,  along  with  the  conversion  to  dR 
(137  MHz)^  is  shown  below. 


N.S.SQ  (/I  = 63°) 

S.  (400  MHz) 
4 

dR  (137  MHz) 

K s 1+ 
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. 0375 

2.  5 

2^  s K « 3+ 
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3.  6 

K > 4 
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. 07  60 

6.  2 

10,  Evans,  .I.V.,  ed.  (1373)  The  Millstone  Hill  Radar  Propagation  Study:  Scientific 
Results,  Lincoln  Lab.  Technical  Report. 
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Figure  A2.  Mean 
Scintillation  Index  for 
Different  Levels  of 
Magnetic  Activity 


This  data  used  by  Evans  was  taken  in  1972  and  1973,  and  was  averaged  over  local 
time  to  include  all  passes.  The  model  dB  values  are  obtained  by  using  the  local 
time  equations  shown  in  Figure  2b  for  each  station  for  three  particular  values. 
A solar  flux  level  of  S^.  = 108  was  used  since  solar  flux  in  1972  and  1973  was  of  this 
magnitude.  The  resulting  averages  over  all  hours  and  all  months  are: 
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Sagamore  Hill  (54°  ^) 
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The  model  dB  values  (converted  to  S^)  are  shown  as  individual  points  in  Figure  A2 
to  facilitate  comparison  witii  Kvans'  data.  The  results  of  Fvans'  data  and  the  model 
dB  values  predicted  are  in  good  agreement  for  both  Narssarssuaq  and  Sagamore 
Hill.  At  Goose  Bay,  the  model  values  are  much  lower  for  the  low  magnetic  activ- 
ity case,  somewhat  lower  for  the  case  K = 2.  5 and  higher  for  the  case  K =5.0. 

P P 

It  is  believed  that  the  position  of  the  trough  is  approximately  the  invariant 
latitude  of  Goose  Bay  under  quiet  magnetic  conditions,  moving  farther  south  as  the 
level  of  magnetic  activity  increases.  The  data  of  fioose  Bay  may  be  affected  bv 
this  phenomenon.  N'either  the  Millstone  Hill  data  nor  the  model  dB  values  have 
taken  into  consideration  the  viewing  geometry.  As  mentioned  earlier,  this  will  be 
part  of  I'hase  2. 
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Appendix  B 


High  Latitude  Models  and  Averaged  vs  Median  Data 

In  Figure  BI  is  shown  the  equation  which  results  in  the  best -fitting  model  for 
N'arssarssuaq  and  Goose  Bay,  for  the  latitude  range  SO"  to  63°A.  The  12  monthly 
plots,  in  local  time,  are  shown  in  Figures  Bla  through  Bll  for  the  cases  of  low  and 
high  solar  flux  (80  and  140)  under  four  selected  magnetic  conditions  (K^  = 0,  5, 

2.  5,  4.  5 and  6.  5),  Under  comparison  with  the  individual  station  plots  exhibited 
earlier,  the  major  cause  of  inaccuracy  in  the  fit  is  seen  to  be  due  to  the  erratic 
high  Kp  model  in  certain  months,  particularly  at  Goose  Bay.  This  is  possibly  the 
result  of  insufficient  data. 

In  order  to  substantiate  the  exclusion  of  Sagamore  Hill  data  in  the  global  model, 
an  additional  equation  was  produced  (Figure  B2)  which  results  in  the  best -fitting 
model  when  data  from  all  three  stations  were  used.  These  equations  are  restricted 
to  the  latitude  range  53°  to  63° A,.  These  monthly  plots  (Figures  B2a  through  B21) 
compare  surprisingly  well  with  those  generated  by  the  two-station  model. 

The  actual  data  plotted  in  Figures  3,  5,  and  10.  upon  which  the  model  is  based, 
are  the  averaged  scintillation  data  in  dB.  It  was  desired  to  compare  this  mean  with 
the  calculated  median  for  the  same  time,  K^,  and  solar  flux  blocks.  Figures  B3a 
B3b,  and  B3c  illustrate  both  the  mean  and  median  dB  values  for  Narssarssuaq, 
Goose  Bay,  and  Sagamore  Hill,  respectively.  This  aata  is  for  the  same  time  per- 
iod (March)  as  described  in  the  earlier  figures.  It  can  be  seen  that  no  substantial 
differences  occur  between  the  averaged  and  median  dB  values  at  any  of  the  stations, 
except  that  as  discussed  under  the  topic  of  distribution  of  SI  readings,  for  low  scin- 
tillation activity  averaged  values  tend  to  be  higher  than  median  values. 
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I.  Narssarssuaq  — Goose  Boy  bosed  . 


-I  3 -2  7CL  + 6.2(I-,I7(  I +2  .6CL)  FD)  [ l + .33(  I -.35 ( I - 8CL)FD)cos(HL  + 2.l-4Kp) 

+ .05cos(2(HL-.9))  + .02cos(3(HL  + 3.5))]  (I-.9CL)  EXA 

. cvA  ,f.2Kp(l  + .l3FD)+.75AsCL(l  + FD)l 
where  EXA  = 2 >■ 


FD  = cos(DA-2.6)  + .2cos(2(DA-30)) 
CL  = cos  ((PF-53.5)  ir/25) 


As  * S^/  100 

HL  = local  time  (hr.) 

PF  = Geomoq.  Lot.  (deg.) 

^ a ^ implied  in  all 
• annual  S diurnal  periodic  terms 
(see  text) 


Figure  Bl.  Mathematical  Model  of  Scintillation  Based 
Upon  Narssarssuaq  and  Cioose  Bay  Data 


Legends  for  data  shown  in  Figures  Bla  through  Bll  and  B2a  through  B21: 

Bla  through  B 11  B2a  through  B21 

Narssarssuaq  Narssarssuaq 

Goose  Bay Goose  Bay 

Sagamore  Hill 
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Plotted  for  March 


Figure  Bid.  Model  Plotted  for  April 


Figure  B le.  Model  Plotted  for  May 


Model  Plotted  for  July 
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H . Narssorssuaq  - Goose  Boy  - Sagamore  Hill  based  . 

SI  = -2.0  + l.2CL  + 6.5(l-.2FD)  [l  + .6(l-.l6FD)(l-.5SL)cos(HL  + 2-.4Kp) 
+ .05cos(2(HL-.9))  + ,02cos(3{HL  + 3.5))]  ( l-CL)  EXA 

[ .2Kp(l  + .2CL)(l-.IXD)  +As(l-.8SL)(I+.45FD)1 
where  EXA  = 2'-  •* 


CL  = cos  ((PF-54  4 ) TT /25  ) CDL  = I + .7SL 
SL=  sin((PF-54.4)  tt/25)  SDL  = I-.3SL 
PD=  cos  (DA+ 16) +.3C0S  (2(DA-30))  ; FD  = CDL»PD  XD  = SDLxPD 

Firjure  B2.  Mathematical  Model  of  Scintillation  Based  I’pon 
Narssarssuaq,  Cloose  Bay,  and  Sagamore  Hill  Data 
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Figure  B2b.  Model  Plotted  for  February 


Figure  B2c.  Model  Plotted  for  March 
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Fif?ure  B3a.  Mean  and  Median  Scintillation  for  yiarssarssuaq  (March) 
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Figure  B3c,  Mean  and  Median  Scintillation  for  Sagamore  Hill 


